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Q^' Abstract 

It is shown that differences observed between the ionization poten- 
tials of the molecular-doped metallic clusters and those corresponding 
to the bare metallic ones can be explained by a two-step approach of 
the classical Liquid Drop Model. This approach takes into account 
' the distinct physical properties of the interface between the molecular 

^ , core and the metallic shell. Also, it is shown that the presence of 

the molecular core may act in the determination of the predominant 
channel of the coulombic fission. 
^ : PACS : 36.40.-C ; 36.40.Wa ; 36.40.Qv 

• 1— I I 

^1 Much of the recent interest in binary clusters has focused on the stabil- 

r-| \ ity of atomic micro-objects consisting of metal atoms surrounding an elec- 

O-i' tronegative impurity There are mainly two interesting 

features exhibited in the behavior of the binary clusters. Most striking is 
^ . the appearance of non-stoichiometric structures, i.e. the octet rule of the 

chemical bonding is apparently broken, in the size range of small clusters (a 
few metal atoms surrounding the electronegative element) 0, ^, ||, ||, 0]. 
Such systems are usually called hypervalent molecules. Nevertheless, by in- 
creasing progressively the number of metal atoms, the transition to metal- 
licity appears as expected. It has been well demonstrated that the pas- 
sage towards the metallic state is accompanied by a structural transition 
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[ly, |Tl], |T2|, |T^, |n|, |1^, |T^, |T^, . Accordingly, in the size range of metal- 
rich clusters, a segregation shows up between a stoichiometric ionically bound 
part and the excess metallic component. 

It has been suggested that the molecular core segregated around the im- 
purity acts to confine the valence electrons into a higher density . This is 
a quantum size effect which results in a perturbation of the electronic den- 
sity of states by a cluster-size dependent amount. Usually, this effect can be 
observed in the ionization potentials [|1^, |1^ |18|, |19|. For example, the ion- 
ization potentials of the lithium monoxide cluster (Li„0) are systematically 
larger (0.1 — 0.15 eV) than the values of the bare lithium clusters (Lin) in 
the size range 2 < n < 70 |l^. Although the difference is much smaller, the 
same trend is also observed for lithium monocarbide clusters (LinC) [|T^]. 
In the context of the quantum size effect, observed differences in the val- 
ues of the size dependent ionization potentials of Li„0 and Lz„C clusters in 
comparison with bare Lz„ clusters have been interpreted and argued for in 
terms of distinct internal contributions (the binding energy and the surface 
potential barrier) to the work function |T9|. A cluster-size dependent work 
function has been derived for doped clusters which is sensitive both to the 
electronegative character of the impurity and to the size of the molecular 
core surrounding the impurity. 

The quantum size effect as we briefly introduced above can be observed, 
e.g. for clusters having spherical shapes, as long as the radial dimensions of 
the molecular and metallic parts are comparable. When the metallic part is 
much greater, the presence of the electronegative impurity is fully screened 
by the valence electron gas and the metallic characteristics of doped metallic 
clusters tend to become similar with those exhibited by bare metallic ones. 



This has been proven in our previous paper [|I9 . 

Nevertheless, the presence of the electronegative element with its sur- 
rounding molecular part may affect the behavior of the binary clusters in 
physical processes which can be explained without involving necessarily their 
quantum aspects. For instance, one can make use of a classical Liquid Drop 
Model (LDM) which, generally, allows an accurate computation of the ion- 
ization potentials and, particularly, offers a suitable description of the frag- 



mentation processes due to the coulombic flssion [^, |21|, |22|, Because 
in the present case one must account for distinct properties of the molecular 
core, the LDM employed in the following is a two-step approach. At the 
flrst instant, one can speculate on the fact that the molecular part segre- 
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gated around the impurity acts itself as an entity by inducing a tension at 
the interface with the metalhc outer shell. Therefore, the molecular core may 
behave differently under shape deformations. This fact certainly changes the 
energy balance between the charged (deformed) state of the doped cluster 
and that corresponding to the bare one. This difference can be observed in 
the ionization potentials [0, |1^ |T8|, |T^. On the other hand, the molecular 



core can equally be seen as a mixing part which gives rise to a configura- 
tional entropy contribution. The presence of the entropic term in the total 
free energy of the system may affect the determination of the predominant 
fission channel for charged (deformed) binary clusters. The behavior of the 
molecular-doped metallic cluster associated with these two particular aspects 
is the focus of this paper. 

The segregation of the molecular inner core changes locally the inter- 
atomic potential and chemical bond lengths. The change of the potential 
energy surface rises restoring forces acting against formation of the curved 
interface between metallic and molecular parts. The fact can be expressed 



in terms of interface tension [El 



(/?) - a,„ (oo) (^1 - ^ j . (1) 

Here, (^^^ is the local curvature, cr^ (oo) is the interface tension correspond- 
ing to the planar interface and 6m is the well-defined microscopic length 
(Tolman's length) p4| depending on the characteristic length of the inter- 
atomic potential inside the molecular part. In the same fashion, the tension 
at the cluster surface can be defined as 

aM (i?o) = (oo) (^1 - . (2) 

All the parameters in above are defined with respect to the properties of the 
metallic part. (-^^ describes now the local curvature of the outer surface of 
the metallic part. 

The radius of the entire cluster Rq can be approximated by 

Rl = R' + {n-p-l)rl , (3) 

where R is the radius of the molecular core, tt, — 1 is the number of the metal 
atoms in the cluster (we considered that one of the initial number n of metal 
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atoms has been replaced by the electronegative impurity), p is the number 
of host metal atoms localized by the electronegative impurity and is the 
inter-electron space (the Wigner-Seitz radii). By specifying the atomic radius 
of the electronegative impurity (ro), the size of the molecular core can also 
be set to 

R' = rl+prl , (4) 

By looking at eq. (1) , one can see that any modification of the local 
curvature of the molecular core changes the interface tension. We can 
ask now the question if the energy balance between the initial and final 
cluster states during the physical processes which involve shape deformations 
maintains the fingerprint of the molecular core. For example, the cluster 
shape undergoes deformations during the ionization process equally due to 
the electrifying phenomenon ||2^ and thermal perturbation Certainly, 



the charged state of a metallic cluster is a deformed one. The deformation 
is presumably over the entire system, which means that the molecular core 
is also affected by this change of the cluster state. It will be shown in the 
following that the information about the molecular core does not cancel when 
considering energy differences between the charged (deformed) and neutral 
states. 

Certainly, the spherical approximation, usually invoked within the jel- 
lium model for clusters |]2T|, becomes less useful for charged clusters whose 



shapes deform due to the balance between electric and surface forces. [20 



Complementary to the jellium model is then the Liquid Drop Model pO| , ^ 



Accordingly, when a spherical cluster of radius i?o, is deformed (conserving 
the volume) toward an ellipsoidal shape, its area changes to 

Am 47ri?^ (l + , (5) 

upon the second-order approximation in powers of the deformation parameter 
(3. In this the simplest deformation, the surface energy is therefore 

Em = 47raM (i?o) {} + ^Z?') • (6) 

Also, in the simplest manner, when the cluster under consideration is elec- 
trified, it has an electrostatic energy 
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where a (r^) is a dimensionless coefficient [p6| . 

As we can see, the surface energy is a minimum when the cluster is 
spherical: The surface energy disfavors deformation. On the other hand, the 
electrostatic energy is a maximum when the cluster is undeformed: It favors 
the deformation. Consequently, the relative strengths of the electric and 
surface forces will determine the shape of a charged cluster which is rather 
deformed. 

As a consequence of the cluster deformation, we can expect that the 
molecular core undergoes deformations, too. The degree of deformation of 
the molecular core depends on the strength of the molecular bonds. In the 
following, we will assume that the change from the sphericity of the molecular 
inner part can be measured by a deformation parameter 7, in the same 
manner as for the entire system 

A„ (7^) - 4rf (1 + ) . (8) 

Consequently, the energy due to the interface between the "molecular" and 
"metallic" parts of the binary clusters becomes 

(7') = 47ra^ (i?) {l + ^7^) . (9) 

At this point one may conclude that the charged state of a cluster is a 
deformed one, the deformation giving rise, in the present case, to changes of 
the surface/interface energies. 

The total energy of a cluster with n atoms and z excess electrons is given 

by 

[ze? 1 



E{n,z) = -zAV + a (r.) ^^l^^^^,/, [l -f')+{n + z) e, + (10) 
+4vraM (i?o) ry/' (l + 

where AV is the outer part of the Coulomb barrier at the cluster surface 
[|T^ and eb is the bulk energy. For a doped cluster, the interface energy term 



given by (9) should be added to the above expression and the size of the 
cluster must be replaced according to eq. (3). The ionization potential is 
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usually defined as a difference in energy between the charged state {z = 



-1) 



WF + a (r. 



neb + 



(11) 



l + >)+ ... , 



where WF = AV — et, and the neutral state of a system {z = 0) 
E{n,0) = neb + A7TaM{Ro)rsn^^^+ ... , 



(12) 



where the deformation of the cluster has been disregarded. (We assume that 
the neutral state of the cluster is spheric.) 

We are now in the position to write the full equation for the ionization 
potentials of molecular-doped metallic clusters. This equation can be written 

as 



IP 



WFb^n + a (r. 



Atxeq [rl + rl{n- 1)] 



1/3 



(l - + (13) 



Stt 



where the distinct physical properties of the molecular inner have been in- 
corporated. The last term comprises corrections due to shape deformations 
of the metallic and molecular parts in the final charged state of the cluster. 
We observe that by replacing in (10) 7 = 0, p = and by the equation 
for ionization potentials for bare metallic clusters is recovered 



IPi 



bare 



WF + a (r. 
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+ -^rl(3'oMni . (14) 
5 



The deformation parameters (3 is inversely proportional with the cluster size 
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^ _ 2 Ar _ 2 Ar 

3 i?o 3 r^ns 



(15) 



where AF is the difference between the semimajor and semiminor axes of 
the ellipsoid. This cancels exactly with the cluster-size dependence of the 
last terms in (14) and makes the surface contribution to depend only on 
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Ar. In the limit of very large clusters this contribution becomes negligible 
(Ar = 0) and, in the asymptotic limit, IPbare —>■ WF. The same applies 
in (13) with regard both to the surface and interface contributions. As 7 
is directly related to /3, the deformation of the molecular part during the 
ionization process being a consequence of the deformation over the entire 
system, the influence of the molecular core will also diminish in the limit of 
large clusters. 

Looking at (13) and (14) , one can see that, apart the surface/interface 
contributions, the other factors which promotes the difference between the 
ionization potentials of the binary, metal-rich and bare metallic clusters are 
the reduction of the valence electron number of the cluster [n — p instead 
of n) and the change of the work function (the usual work function WF 
is replaced in (13) by a work function for binary clusters, that is WF^in) 
[|r^. This change of the work function for doped clusters is a quantum 
size effect due to the perturbation of the electron density. This difference 
disappears with increasing cluster size, the impurity being fully screened by 
the remaining free electrons for a cluster with about 100 atoms, as it has 



been proven in our previous work |]T9 



In the following, we disregard the distinct internal contributions to the 
work function of doped clusters arising from the quantum size effect {WF = 
WFbin), and focus on the surface/interface contributions promoted by dis- 
tortions of the cluster shape. We derive the ratio between IP and IPbare , 
which is 



IP 



IP 



bare 




(16) 



We make a numerical test for the above result and compare with the ex- 
perimental observation as reported in Ref. [0. For example, in the size 
domain of of large clusters (n > 6), the experimental IP's of oxygen doped 
lithium clusters are systematically larger (about 4%) than the values of the 
bare lithium clusters. By looking to the above equation (16), we see that for 
n ^ 1, this can can be approximated by 

J^^l + t^fP-Y' (17) 

IP,are fi^aM\n) ^ ' 
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which leads to 

We set p = 2 and n = 20 (this is supposed to be the case of (Lz20) Liig 
cluster) and observe that the difference IP — IPtare of about 4%, as reported 
in Ref. |TP[, can be obtained if either the tension is higher at the surface 
than at the interface (ctm — 5am) at 7 = /3, or, equivalently, ^ < 0.04 at 
(Tm — cTm- The latter means that the molecular core is less affected by the 
distortion of the entire cluster. The molecular core opposes to the external 
action exerted by the metallic outer shell. Indeed, both effects can concur to 
set the above difference. 

At small cluster sizes, when the charged clusters are rather deformed, the 
interface contribution in (9) can play an important role. This has a material 
dependence by a^n (R)- There is a sizeable body of work on aspects of the 
statistical mechanical theory of interfacial phenomena [^]. Anyway, it is 
difficult to say what results tell us about the physical properties of a real 
interface. The concept of a bare surface/interface tension is controversial, 
since correlation effects for the restoring force contribute also to the specific 
thermodynamical potential. Looking at (1) , we can see that the curvature 
affects the surface tension. The effect comes to light at full extent for micro- 
scopic drops whose radii approach the range of intermolecular forces. The 



Tolman length S entering (1) has negative values |24 

6 ^ -0.39C , 

which depends on the characteristic lengths of the intermolecular potentials, 
^. Obviously, we can imagine that different properties of the binary clusters 
can be explained by taking the particular chemical structure of the molecular 
part into account. Some geometrical structures of the molecular part can 
fetter the ellipsoidal deformation of the entire system or, by contrast, others 
can favor it. So that, the geometrical structure of the molecular part can 
play a role in moderating the distortion tendency of the charged clusters. 

Another aspect related to the presence of the molecular inner core is 
the behavior of the binary clusters during the fragmentation processes by 
coulombic fission. It was shown that, generally, the binary fission occurs 
predominantly with asymmetrical character |2^. The discrete nature of the 
cluster constituents become important in this case, since the same fragments 
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may consist of different atoms. In order to see how the molecular core can 
affect the determination of the predominant fission channel, we describe the 
fission process 



matter of minimizing the free energy PH EUI 

F = T, + C -keTS , 



where 



E 
C 



aj^q'^/^ + af^_^ (n - q) 

^2 



(19) 



(20) 



CnU 



-1/3 



-1/3 



-1/3- 



are the surface and coulombic energies, a/- {i = q,n — q) and a„ stand for 
the surface energies per atom for the fission fragments and parent sample, c 
is the coulombic energy per atom and the subscripts have the same meaning 
as above. The entropy 5" is given by 



ks {n-q)\q\ 



(21) 



-n 



n n 



'^^Infl-^ 



n 



n 



and shows that the same fragments may consist of different atoms. We 
notice that in the extreme case when one of the fission product, let say q, 
is just the molecular core, the surface (a/,) and coulombic (cf^^ energies, 
which control the fission process, must be properly replaced in (20). The 
surface term corresponding to the fragment q can be derived starting from 
(1). The estimation of the coulomb energy of the fragment q must account 
for the appropriate character of the chemical bond in the molecular core. To 
compare with the trivial case of bare metallic clusters, we make the terms 
Qf. and ttp equal each other (a/. = Op = a), and replace c/^ = | c/^-, = i 
and Cp = Y [iH • We observe that the predominant fission channel for binary 
clusters is affected by the physical properties of the molecular core. We can 
also see that if the molecular core is more stable than the rest of the metallic 
cluster then, the size of the molecular core acts as a lower bound for the 
characteristic mass of the fission product {q > p + 1) |T^. 
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Finally, the present study on binary, metal-rich clusters based on a two- 
step LDM approach allowed us to point out their different behavior in com- 
parison with the corresponding bare metallic species. This behavior has a 
general feature for metal-rich clusters, in agreement with the segregation ev- 
idence of the molecular part: The molecular core acts itself as an entity by 
inducing a tension at the interface with the metallic outer shell and as a mix- 
ing part by giving rise to a configurational entropy contribution. The former 
depends on the chemical character of the binding and its energetic contribu- 
tion is in direct proportion with the interface area. The latter comes to light 
at full extent for rather high temperatures. The effect of the former leads 
to an energetic balance in the final (charged) state of the cluster wherefrom 
the difference between the ionization potentials of the doped clusters and 
corresponding pure clusters shows up. The latter affects the determination 
of the predominant fission channel of the coulombic fission. The size of the 
molecular core can act as a lower bound for the characteristic mass of the 
fission product. 

Concluding, one can say that, although the metallic characteristics def- 
initely predominate in the range of large cluster sizes, the electronegative 
impurity still affects (not dramatically, indeed) their thermodynamical prop- 
erties. Thus, the molecular core of the binary, metal-rich cluster has a re- 
sponse to external actions, this resulting in a change of the total energy of 
the system. 
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